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Abstract. We present an analysis of structural features 
of all known galaxies with optical polar rings. We find a 
clear dichotomy for objects of this peculiar class. Bulge- 
dominated SO galaxies possess only short narrow rings, 
while disk-dominated objects always have wide extended 
polar rings. We try by gas dynamical simulations to ex- 
plain such a segregation by dependence of the ring-forming 
process on different galaxy potentials. It is found that the 
total mass captured into the ring during an encounter of 
a host-ring system with a gas-rich spiral galaxy of com- 
parable mass exceeds 10 9 M Q (or about 10% of all gas in 
the donor galaxy), which is of the order of that found by 
observation. The process of gas to gather into a steady- 
state ring takes approximately (7 — 9) x 10 8 years. This 
time is somewhat shorter for rings forming around bulge- 
dominated galaxies. We also present observational argu- 
ments for SO galaxies with extended rings to be similar 
to late-type spirals by their photometric properties, while 
numerical modelling of the extended ring formation sug- 
gests that these galaxies must possess massive dark halos 
as well. In this case, the sizes of the modelled rings turn 
out large enough (up to 30 kpc in diameter) , and the time 
scale for ring formation is prolonged up to several Gyrs. 

Key words: galaxies: general, interactions, photometry, 
kinematics and dynamics, peculiar, structure - dark mat- 
ter 



1. Introduction 

Recent observations (and especially with HST) demon- 
strate that mutual interactions and mergers between 
galaxies at early stages of evolution of the Universe were 
probably among the main processes leading to the ob- 
served properties of galaxies (e.g., Keel 1996). Even at 
the present epoch, at least 5-10% of galaxies are members 
of interacting systems. Many other galaxies keep signs 
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in their structure of past interactions and mergings (for 
example, elliptical and SO galaxies with inclined gaseous 
disks, galaxies with faint shells and ripples, galaxies with 
kinematically decoupled nuclei, etc.). 

Polar-ring galaxies (PRGs), consisting of large-scale 
rings of stars, gas and dust orbiting around major axes 
of early- type galaxies, may be considered as extreme sam- 
ples of possible interaction relics. Indeed, in the case of 
PRGs, the remnants of merged galaxies are not mixed in 
one smooth object but stay separately in a quasi-steady 
state for a long time. PRGs are very rare objects. For 
example, the Polar Ring Catalogue by Whitmore et al. 
(1990) (PRC) lists only 6 classic kinematically-confirmed 
polar-ring galaxies. 

The unique geometry of PRGs attracts the attention 
of astronomers trying to test the 3D shape of galactic po- 
tentials and to study the stability of inclined rings and 
disks (see recent review articles by Tohline 1990, Sackett 
1991, Sparke 1991, Combes 1994, Cox & Sparke 1996). 
Such an important question as the origin of these peculiar 
objects it is still not adequately investigated. It is usu- 
ally suggested that collapse of a single protogalactic cloud 
cannot create an object with two nearly-orthogonal large- 
scale systems (but see Curir & Diaferio 1994), and so some 
secondary event must occur in the history of PRGs. 

Summarizing possible scenarios of polar-ring forma- 
tion, one can enumerate the following: the accretion of 
matter from a nearby system or the capture and merging 
of a gas-rich companion; the delayed inflow of a primor- 
dial intergalactic cloud; the accretion of matter from the 
outskirts of the host galaxy itself; the polar-ring formation 
from the return of the tidal material during the merging of 
two gas-rich spirals (Toomre 1977, Shane 1980, Schweizer 
et al. 1983, Sackett 1991, Sparke 1991, Hibbard & Mihos 
1995). 

Probably, all the above mechanisms can create ring- 
like structures around early-type galaxies. To our mind, 
the most straightforward scenario is the first one. Recent 
observations of several binary interacting systems clearly 
demonstrate such rings in the making (for instance, NGC 
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7464/65 - Li & Seaquist 1994, NGC 3808A,B and NGC 
6285/86 - Reshetnikov et al. 1996). 

Another unclear question is the nature of central ob- 
jects in PRGs and a possible correlation of host galaxy 
properties with characteristics of a surrounding polar ring. 

In this paper, we present SPH simulations of polar 
ring formation around target galaxies of different struc- 
tures due to gas accretion during the encounter with a 
comparable-mass spiral galaxy. In our simulations, we try 
to answer the following main questions: Does this mech- 
anism work? What determines the size of the resulting 
ring, and what is its spatial structure? On what timcscalc 
does the ring form? What is the mass fraction of the gas 
captured into a ring? 

The paper is organized as follows: in Section 2, we ex- 
amine observational properties of all known kinematically 
confirmed PRGs and formulate some observational con- 
straints on numerical simulations; in Section 3, we discuss 
previous attempts to model the PRGs formation, describe 
our modelling technique and results of simulations; and 
finally we give our conclusionsin Section 4. 

Throughout the paper, all distance-dependent quanti- 
ties are calculated using Ho = 75 km/s/Mpc. 

2. General characteristics of polar-ring galaxies 

As a definition of a polar-ring galaxy, we will use the def- 
inition of Category A objects in the PRC: spectroscopic 
evidence must exist for two nearly-perpendicular kinemat- 
ical subsystems; centers of the two components must be 
aligned, and both subsystems must have similar systemic 
velocities; the ring must be comparable in size to the host 
galaxy, must be luminous and nearly planar. This defini- 
tion allows to separate dust-lane ellipticals, galaxies with 
inclined HI rings etc. from PRGs. Using this rigorous def- 
inition, one can now consider only three additional galax- 
ies to 6 the classic PRGs listed in the PRC: AM 2020-504 
(Whitmore & Schweizer 1987, Arnaboldi et al. 1993), IC 
1689 (Reshetnikov et al. 1995, Hagen-Thorn & Reshet- 
nikov 1997) and NGC 5122 (Cox et al. 1995). (We do not 
consider ESO 603-G21 here due to the puzzling kinematics 
of the central galaxy (Arnaboldi et al. 1995).) 

An examination of the optical images of PRGs (e.g. 
in the PRC) allows one to divide them into two groups 
(Whitmore 1991): galaxies with extended disk-like rings 
with the central region cut out and galaxies with relatively 
narrow rings, not extended in radius. This division is quite 
distinct since the first group of galaxies - A0136-0801 (A- 
1), UGC 7576 (A-4), NGC 4650A (A-5), UGC 9796 (A-6), 
and NGC 5122 (B-16) - possess optical rings extended out 
to 2-3 diameters of the central galaxies, while the second 
group - ESO 415-G26 (A-2), NGC 2685 (A-3), IC 1689 (B- 
3), and AM 2020-504 (B-19) - demonstrate optical rings 
with size not exceeding the diameter of the host galaxy. 

In Table 1, we generalize the main observational char- 
acteristics of the two groups of PRGs. (Note that, due to 



the absence of optical data about NGC 5122, we did not 
consider this galaxy in the table.) In the case of incom- 
plete data or large scatter of characteristics, we give in 
the table only limits or indicate the range of parameter 
changes. Absolute luminosities and colors in the table are 
corrected for Galactic absorption. 
We discuss Table 1 in detail. 



2.1. Central galaxies 

The problem of the bulge-to-disk ratio determination for 
the PRGs is not quite simple. Reshetnikov et al. (1994) 
have noticed that bulge effective parameters (/i e and R e ) 
of PRGs with wide extended rings (our first group) lie 
in the plane of effective parameters below the mean rela- 
tionship for normal galaxies (this means that at the same 
value of /i e , bulges of PRGs have smaller radii in com- 
parison with bulges of normal galaxies) - see Figure 7 in 
the cited paper. But bulge characteristics of the second 
group of galaxies are located exactly along the standard 
relationship. Such a dichotomy suggests two explanations: 
the unusual compactness of the first group of PRG bulges 
or the underestimation of their sizes and luminosities due 
to the projection of gas and dust rich extended rings on 
central regions of the galaxies (Reshetnikov et al. 1994). 

A comparison of photometric cuts of the galaxies in 
different color bands allows one to solve this problem. In 
Figure |l|, we present surface brightness distributions for 
four PRGs with extended rings along the major axes of 
their central galaxies. We compare two profiles for each 
galaxy: one in the B passband and other in the red filter 
(R, i, or K). As one can see from figure, absorption in the 
rings barely changes the shape of the profiles. Therefore, 
PRGs with wide extended rings in fact possess unusu- 
ally compact and faint bulges in comparison with normal 
early- type galaxies. As is evident from Figure [l], the expo- 
nential disk dominates in the photometric structure of all 
these galaxies. The ratio of the total bulge luminosity to 
the luminosity of the exponential disk in the first group of 
PRGs is about 0.1 (we denoted this as << 1 in Table 1). 

The photometric structure of PRGs with short rings 
(the second group of galaxies) looks usual for early-type 
galaxies. As was mentioned earlier, their characteristics in 
the plane of effective parameters follow the mean relation 
for normal galaxies. According to original papers, bulge- 
to-disk ratios in the B band for these galaxies are: ESO 
415-G26 - 0.52 (Whitmore et al. 1987), NGC 2685 - 0.9 
(Makarov et al. 1989), IC 1689 - 1.9 (Reshetnikov et al. 
1995). AM 2020-504 is an elliptical galaxy (Arnaboldi et 
al. 1993). Therefore, the second group PRGs are normal 
bulge-dominated galaxies by their photometric structure 
(conditionally, « 1 in the table). 
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Parameter Extended rings Short rings 
Central galaxy: 

Bulge-to-disk ratio << 1 ~1 

Absolute luminosity in the B band -18.6 ± 0.5(a) -19.2 ± 0.3(a) 

B-V +0.84 ± 0.05 +0.91 ± 0.03 

Colour gradient Yes Yes 

a (km/s) 80 ± 20 100 + 270 

V m ax/a 1.7 ± 0.3 0.45 + 1.4 

Decoupled nuclei ? Yes 

Ring: 

Diameter « (2 - 3)D 25 < D 25 

26 ± 6 kpc 9 + 4 kpc 

Absolute luminosity in the B band —17.5 ± 0.6 > —15.4 

B-V -0.1 to +0.7 +0.2 to +0.64 

Colour gradient Yes ? 

HI mass (1O 9 M ) 6 + 2 < 9 

Ring inclination 7° + 26° 1° + 16° 



Fig. 1. Photometric profiles of PRGs with extended rings along major axes of the central galaxies. The data for A0136-0801 
are from Schweizer et al. (1983) and PRC; NGC 4650A - Whitmore et al. (1987) and Combes & Arnaboldi (1996); UGC 7576 
and UGC 9796 — Reshetnikov et al. (1994). The dashed lines represent exponential fits of the observed distributions. 



2.2. Polar rings 



Rings sizes in the two groups of galaxies are different, 
as well as their absolute luminosities (D25 in Table 1 de- 
notes a diameter of the central galaxy measured at the 
surface brightness level [Ib =25). Both groups of rings 
show a large scatter of optical colors although, as it was 
shown by Reshetnikov et al. (1995), there is a general 
trend with a blue surface brightness; bluer rings have, on 
average, higher surface brightnesses. Extended rings also 
show large-scale color gradients, that is, they become bluer 
at larger radii (see, for instance, Figures 4,6 in Reshet- 
nikov et al. 1994 and Arnaboldi et al. 1995). Both group 
of rings contain a large amount of neutral hydrogen (we as- 
sume in Table 1 that all detected HI belongs to the rings), 
with larger scatter of HI mass in PRGs with narrow rings 
(Schechter et al. 1984, van Gorkom et al. 1987, Richter 
et al. 1994). The last line of Table 1 shows the range of 
ring inclinations (angular distance between the ring and 
the perpendicular to the central galaxy plane) according 
to Whitmore (1991). Both groups of PRGs demonstrate 
rings to be very close to perpendicular, with a somewhat 
larger deviation for extended rings. It should also be noted 
that global characteristics of extended rings resemble the 
disks of spiral galaxies (this was previously pointed out 



by Reshetnikov et al. 1994 and Reshetnikov & Combes 
1994a,b). We will discuss this analogy in Section 4. 

Summarizing the above analysis, one can conclude that 
there is a correlation between general properties of optical 
polar rings and characteristics of host galaxies. Extended, 
disk-like rings exist preferably around galaxies which have 
global photometric structure (and probably mass distri- 
bution) similar to late-type spiral galaxies. Host galaxies 
of PRGs with extended rings demonstrate a remarkable 
similarity — their characteristics fall in a relatively narrow 
range (see Table 1). From the other side, bulge-dominated 
galaxies have tighter and narrow rings. Host galaxies and 
rings of PRGs with non-extended rings demonstrate sig- 
nificantly larger scatter of intrinsic properties. 

Although the statistics of PRGs are not great (8 - 
9 objects only), the mentioned tendency is quite dis- 
tinct — there are no bulge-dominated galaxies with ex- 
tended luminous polar rings. It is significant since the 
known PRGs were selected on the base of optical mor- 
phology only, without any attention to the structure of 
the central galaxies. One can note also that the PRG di- 
chotomy is quite analogous to a recently-found difference 
of ionized gas distribution in elliptical and SO galaxies. 
According to Macchetto et al. (1996), more than half of 
elliptical galaxies with detected ionized gas have their gas 
concentrated in small (<4 kpc) nuclear disks, while most 
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SO galaxies demonstrate more extended (up to 18 kpc) 
distribution of the gas. 

A similar tendency was mentioned for the first time by 
Whitmore (1991), who found that only rapidly- rotating SO 
disks have extended luminous polar rings, while dust- lane 
ellipticals rarely show an extended luminous component. 
He suggested that a more flattened potential of the SO 
galaxy is able to stabilize the ring at greater radii than in 
an elliptical galaxy. In the next Section, we will demon- 
strate that the observed dichotomy could be caused by 
another reason. 



3. Numerical simulations 

There are only a few investigations that illustrate several 
proposed scenarios of polar ring formation. This lack of 
theoretical study devoted to the question of ring origin 
is connected with the necessity of the usage of compli- 
cated 3D gas dynamics models, since polar rings contain 
a high amount of neutral gas (up to a few 10 9 M ), as seen 
from Table 1. One of the most fit tools to construct such 
models is a smoothed particle hydrodynamics algorithm 
- SPH, - originally proposed by Lucy (1977) and Gin- 
gold & Monaghan (1977) and considerably developed in 
the late eighties (e.g., Hernquist & Katz 1989). It is fully 
three-dimensional, flexible regarding the geometry of the 
gas distribution and consequently well-suited to modelling 
both the donor and recipient galaxies and the empty ex- 
tended space between them. Applied to a range of polar- 
ring investigation problems, this method has already lead 
to many important results. 

The main result of the ring stability investigations is 
that the time for gas in an axisymmetric or triaxial po- 
tential to settle into a steady state is small compared to 
the age of the Universe (see Christodoulou et al. 1992 and 
references therein) . So there is sufficient time for accreted 
gas to form polar rings. However, these studies did not 
take into account time-dependent effects and realistic ini- 
tial conditions. Some efforts to include these effects have 
been made by Rix & Katz (1991) and Weil & Hernquist 
(1993). 

Rix & Katz (1991) have treated the polar-ring for- 
mation process as a gradual smearing of a small gaseous 
blob captured by a galaxy of spherically-symmetric struc- 
ture on a distant circular orbit. Naturally, the diameter 
of the forming ring (up to 80 kpc) was quite completly 
dictated by the initial position of a gaseous satellite (the 
orbit radius). The approach of Weil & Hernquist (1993) 
was more self-consistent. They have considered a parabolic 
encounter of a low-mass gas-rich companion with a more 
massive galaxy and consequent merging of the former. The 
gaseous component of the satellite was settling into the 
ring after the total disruption of the satellite in the vicinity 
of the target galaxy. The forming ring turned out rather 
small (about 6-7 kpc in diameter) and its size was de- 



termined in the end by a hand-introduced description of 
details of the companion disruption. 

The ability of the accreting material to dissipate en- 
ergy is the crucial factor for ring formation. Test (nonin- 
teracting) particles (which imitate molecular clouds from 
intcrgalactic space or clouds belonging to a companion 
galaxy) lose energy due to ram-pressure in an extended 
halo of a spiral galaxy and form a gaseous ring rotating 
around the disk of the spiral (Sofue & Wakamatsu 1993, 
Sofue 1994). 

3.1. Modelling of polar-ring characteristics during the 
ring-forming process 

The main feature of the present investigation is a descrip- 
tion of the full history of the gas stripping of the spiral 
galaxy outskirts and its consequent capture by a satellite 
during a parabolic encounter. Keeping in mind the recent 
observations of several interacting pairs of galaxies of com- 
parable luminosities which demonstrate evidently ring- like 
structures in the making (see Introduction), we have con- 
sidered a distant encounter of equal-mass systems. In the 
case of a distant encounter the effect of self-gravity may 
be ignored. The treatment of the gas hydrodynamics is 
described below. 

We did not explore in a comprehensive manner the or- 
bital parameter space and choose only one set of impact 
parameters for which the modelling encounter of an SO 
galaxy with a gas-rich system of comparable mass unam- 
biguously results in the polar-ring formation around the 
former. Motivated by the above observational data anal- 
ysis, we have undertaken a numerical investigation of the 
ring-forming process for galaxies with different structures. 

3.1.1. Method 

Our simulations are based on a rather standard variant 
of the SPH code, which is the same as that in Sotnikova 
(1996). For simplicity, we adopt the smoothing length h 
(analogue of the particle size) to be independent of the 
local gas density. The smoothed values of hydrodynamical 
parameters are estimated by using a grid with cell length 
equal to 2h. We assume an isothermal equation of state, 
leaving thermal effects completely unexplored. The gas 
is always at a temperature 10 4 K, and the corresponding 
equation of state is 

P = c 2 Pl 

where P and p arc the gas pressure and density, and 
c = const is the isothermal speed of sound. For 
T = 10 4 K , c w 9.1km/s. 

The motion of N = 10 000 particles, which represent 
elements of the continous gaseous medium, is traced by 
means of equations of momentum conservation: 

dr t 

-dt =Vi > W 
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dt 
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where and Vi are the spatial coordinate and velocity of 
the z-th particle, and $ is the gravitational potential. The 
term a^ lsc describes the acceleration due to viscosity qc 



1 

Pi 



where pi — p(ri). We adopt one of the standard forms 
for the artificial viscosity qi which serves to represent the 
sudden deceleration of the gas motion when strong shocks 
arise in the gas (e.g., Hernquist & Katz 1989): 



= (ach Pi \V ■ Vi\ + (3h 2 Pi \V • ^| 2 , V • v { < ^ 
[0, otherwise . " 



In the expression (||), the viscosity depends on the diver- 
gence of the velocity field. The first term is analogous to a 
bulk viscosity, the second, indroduced to prevent particle 
interpenetration at high Mach number, is of the Neumann- 
Richtmyer type. Parameters a and j3 are free parameters. 
According to Hernquist & Katz (1989), one can satisfac- 
torily reproduce the change of the density and pressure 
across the shock front if the values of a and (3 are equal 
to 0.5 and 1.0, respectively. We used just these numbers 
in our simulations. 

The procedure of a smoothed-value estimation of the 
gas density is quite usual, as well as the transition from hy- 
drodynamical equations (1) and (2) to SPH-equations (for 
a detailed description of the adopted numerical scheme, 
see Sotnikova 1996). 

3.1.2. Model 
Donor galaxy 

To minimize the number of free parameters, we chose 
a very simple model for the donor galaxy - its potential 
was taken as that of a softened point mass 



*o(r) 



GM 



(r 2 +a2)V2 ' 



(4) 



where G is the gravitational constant, Mq the mass, and 
ao the softening scale length of the potential. 

N = 10 000 particles are used to represent the gaseous 
medium. They are initially gathered in a thin disk with 
a 1/r distribution from the center of the donor galaxy up 
to an outward radius Rg and are placed in dynamically- 
cold circular orbits. The total gas mass is 0.2 M . The 
gaseous particle size h, which determines the numerical 
spatial resolution, is 300 pc. 

Throughout, we employ the following system of units: 
the gravitational constant G = 1, the mass of the donor 
galaxy Mq = 1, the outer radius of the initial gas distri- 
bution in the donor galaxy Rq = 1. In terms of physical 
units, a unit mass can be defined as 10 n M©, and a unit 



distance corresponds to 15 kpc. Combined together, they 
give a unit time to be 86.6 Myrs, and a unit velocity to 
169.3 km s _1 . 

Accreting galaxy 

We have considered the two-component model of a 
galaxy, consisting of a bulge and a disk. The light dis- 
tribution of the bulge is well-provided by the model used, 
for example, in Weil & Hernquist (1993). For this model, 
the potential is 



$ b (r) = - 



GM h 
r + a h 



(5) 



where Mb is the bulge mass, and ah the scale length of 
the potential. 

The gravitational potential of the disk component of 
the galaxy is approximated by a Miyamoto & Nagai (1975) 
potential 



GM A 



Jy 2 + z 2 + (b d + y/x 2 + aff 



(6) 



where Md is the disk mass, and ad and bd are scale radii. 
This model was favored by the simplicity of its practical 
usage. In our simulations, a^/bd = 0.2. 

The ratio Mb /Mi determines the range of models from 
bulge-dominated systems (large values of the ratio) to 
disk-dominated galaxies (small values of the ratio) . Four 
models for an accreting galaxy have been considered, with 
values of M h /M d equal to: 2.0, 1.0, 0.5, 0.2. The total mass 
of the galaxy is taken the same for all models and equal 
to that of the donor galaxy: M b + M d = 1 = 10 n M Q . 

As shown in Section 2, main galaxies of PRCs with 
extended rings are similar to spiral galaxies in global 
photometric structure (and probably mass distribution). 
According to de Jong (1996), the bulge effective surface 
brightness shows the best correlation with morphological 
type of spiral galaxies. Therefore, the total bulge lumi- 
nosity (mass) also correlates well. Scaling radii show a 
large scatter for all types of spiral galaxies. Thus, we fixed 
the scale lengths of galaxy components, leaving only the 
mass of the bulge to be variable. Scale lengths of the bulge 
Ob and disk bd potentials are 1 kpc and 5 kpc, respec- 
tively. For these values, the half-mass radius of the bulge 
is (1 + V2)ab — 2.4 kpc, and the maximum of the ro- 
tation curve of the disk component falls on the radius 
V%(bd + ad) ^ 8.4 kpc. 

Orbit 

Before modelling a ring-forming encounter, we first 
solved the two-body problem for the donor and ring-host 
galaxies. The numerical procedure employed is quite sim- 
ilar to that used in Weil & Hernquist (1993). 

Two interacting galaxies are initially separated by a 
rather large distance (rj n ; = 5 = 75 kpc), so that tidal 
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effects be negligible. The initial velocities are taken as for 
parabolic encounter. The primary galaxy passes in a zero- 
inclination (in the plane of the gaseous disk of the donor 
galaxy - xy plane), prograde orbit (that is, the orbital 
angular momentum is parallel to the spin of the donor 
galaxy) with a pericenter distance r por initially set as 1.6. 
The polar axis of the ring-host galaxy disk component 
(parallel to the x-axis) lies in the orbital plane, so the orbit 
is polar regarding the ring-host galaxy. The calculated or- 
bit somewhat differs from parabolic (due to non-Keplerian 
potentials of galaxies), but its form is nearly the same for 
all models of a ring-host galaxy. 

3.1.3. Results 

During the encounter, the primary galaxy strips the out- 
skirts of the donor object and a ring, rotating in the di- 
rection of the orbital motion, eventually forms around the 
galaxy in the encounter plane. As the equatorial plane of 
the disk component of SO is taken to be perpendicular 
to that of the orbital motion, the ring is polar. The to- 
tal amount of accreted gas is about 10% of all gas in the 
donor galaxy (or about 2 x 10 9 M Q ) and there is not any 
significant difference in the stripped mass for all consid- 
ered models of the ring-host galaxy. Some amount of the 
gas fell into the central part of galaxy. The mass of this 
gas is estimated as 2 x 1O 8 M within 1 kpc from the cen- 
ter. The timescale of the ring formation is a few 10 8 years. 
This time is somewhat shorter for rings forming around 
bulge-dominated galaxies and reaches up to ~ 9x 10 8 years 
for a disk-dominated model. The interacting galaxies are 
separated by a sufficiently large distance (more than 120 
kpc) by the time of the ring steady-state settling, so that 
there are no any direct evidences of interaction between 
the galaxies and the ring-host system as an isolated object. 

The simulations have revealed an interesting feature 
which appears during the ring-forming process in different 
galaxian potentials. As a steady state sets in, ring sizes 
begin to diverge. The results for four runs are presented 
in Figure 0. 

Figure |2j gives a final stage of polar-ring formation for 
all considered models. This is the orbital plane projection 
of the gaseous rings. The marked difference in the ring 
structure is the difference in ring size, which rises when 
the bulge mass decreases and ranges from about 7 kpc 
in diameter for bulge-dominated systems (see Table 1) to 
approximately 13 kpc for disk-dominated objects. 

The ring radius is finally determined by the angular 
momentum of donor galaxy gaseous particles forming the 
ring with respect to the galaxy at the moment of the peri- 
center passage. As the impact parameters are nearly the 
same for all runs, the value of the total angular momentum 
of ring-forming particles is almost the same for all models. 
But one can see from Figure || that the positions of par- 
ticles with the same angular momentum are different for 
galaxies with different structures. Particles are closer to 



the center of galaxy if there is a more concentrated mass 
distribution - the bulge-dominated model. This implies, 
in particular, that under the same conditions rings form- 
ing around elliptical galaxies with a strong concentration 
of mass to the center (see, for instance, the curves marked 
as (* — * — *) and (o — o — o) in Figure]^) should be less 
extended (on the average) - rather internal - than those 
around galaxies with a more gently sloping density profile 
(this fact was first mentioned in Sotnikova 1996). 



Fig. 3. The angular momentum of a test particle on a cir- 
cular orbit versus the distance from the center of a galaxy 
in the equatorial plane of the disk for 4 bulge-disk mod- 
els. The short, medium and long dash lines conform to 
Mb/Md = 2.0, 1.0, and 0.5, correspondingly; the solid line gives 
Vrotr — r relation for Mb /Mi =0.2. Two additional lines are 
shown as the bench-marks: (* — * — *)- for diskless model, 
and (o — o — o) - diskless model with more concenrated bulge 
(ab = 0.5 kpc). The horizontal line marks the value of the 
mean angular momentum per ring particle for 4 runs. 

In Figure |4|, the surface density profiles along the ma- 
jor axes of the rings are shown. One can see one more 
interesting feature of simulated rings: rings farther from 
the center of the galaxy are more extended. 

These results give the post factum justification of our 
impact parameter choice. Indeed, a closer passage re- 
sults in less extended rings. The angular momentum of a 
gaseous particle at the moment of the closest approach of 
galaxies with respect the ring-host system is significantly 
affected by the term proportional to ~ v peT r veT , where 
v peI is the ring-host galaxy velocity relative to the donor 
galaxy. For parabolic encounters, it leads to ~ ^/r pcr . As 
the ring size monotonically decreases, the angular momen- 
tum decreases; a closer encounter gives a ring lying rather 
inside the galaxy. It is known that SO galaxies possess some 
original gas, so the interaction between two gaseous sys- 
tems would result in ring destruction. The further decreas- 
ing of pericenter distance leads to merging of the galaxies 
and total change of their structures. As to the more dis- 
tant encounters, they do not allow capture of substantial 
amounts of gas, due to the sharp decrease of disturbing- 
forces. We have repeated our runs with r pcr = 2.4 = 36 
kpc and have obtained that the mass of the resulting 
very diluted rings did not exceed 2 x 10 8 Af Q . Hyperbolic 
galaxy encounters are not efficient for ring formation, as 
the time of strong galaxy interaction is short. Let us note 
also that we discuss the properties of classic PRGs with 
optical star-forming rings. So we discuss the origin of rel- 
atively dense rings - with gas density large enough for 
star formation. This gives an additional restriction on the 
impact parameter, since the matter captured during the 
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Fig. 2. Time 10 (after the pericenter passage moment) frame for 4 runs. The values of the bulge to disk mass ratio of the 
ring-host galaxy are shown in the upper right-hand corner of all frames. All frames display the orbital plane (xy-plane) spatial 
projection and measure 30 kpc per edge. The equatorial plane of the ring-host galaxy disk component (parallel to j/z-plane) is 
perpendicular to that of the orbital motion. 



Fig. 4. Surface density profiles of rings. The short, medium and long dash lines conform to Mb/Md = 2.0, 1.0, and 0.5, corre- 
spondingly; the solid line - Mh/Md — 0.2. 



distant encounter will spread along a more extended orbit 
and will have lower mean density. Thus, galaxy encoun- 
ters favourable to extended optical ring formation are ex- 
tremely rare events. Such interactions must be between 
galaxies of specific types and within restricted geometry 
(polar encounters and a very narrow range of impact pa- 
rameter). 

We found in our simulations that, with the same 
impact parameter, the rings forming around bulge- 
dominated galaxies are less extended (by factor two) than 
the rings forming around disk-dominated galaxies. There- 
fore, having the same accretion history in the samples of 
elliptical and disk-type galaxies, we will have on average 
more extended rings around disk-dominated galaxies. This 
inference may explain in a quite natural manner the ab- 
sence of extended polar rings around elliptical galaxies. 
It can also explain a difference in morphology of ionized 
gas in elliptical and SO galaxies (Macchetto et al. 1996), 
assuming external origin of gas in these galaxies. 

One can remark, however, that the obtained segre- 
gation of ring sizes (see Figure ||) is not as clear as the 
observed dichotomy (see Table 1). The simple bulge-disk 
model giving internal rings for bulge-dominated objects 
fails to explain the existence of very extended (up to 30 
kpc in diameter) polar rings around disk-dominated galax- 
ies. Probably, there exists one more physical factor leading 
to the strong differences of the PRG structural properties. 

3. 2. Exploration of dark halo inclusion effect 

As it was found from the observational data analysis, there 
exists for PRCs with extended rings a remarkable struc- 
tural resemblence of host galaxies to late- type spiral galax- 
ies. It is known that dynamical properties of late-type 
spiral galaxies are determined to a considerable extent 
by invisible massive halos (e.g., Freeman 1992). Hence, 
we can suppose that PRGs with extended rings possess a 
third global component (besides the bulge and the disk) 
- a dark massive halo. The existence of dark massive ha- 
los also follows from the analysis of ring kinematics (e.g., 
Schweizer et al. 1983, Reshetnikov & Combes 1994a). To 



account for the rotational velocities in the rings of UGC 
7576 and UGC 9796 inside radii of 17 kpc and 21.4 kpc, re- 
spectively, one requires amounts of dark mass reaching 1.6 
and 3 times the luminous masses (Reshetnikov & Combes 
1994a). 

What are the possible consequences of taking into ac- 
count this structural component in the context of our in- 
vestigation? As the gravitating mass of the primary galaxy 
increases, the orbital velocity of this galaxy relative to the 
donor object rises also. Then angular momentum argu- 
ments for a galaxy with such a structure lead to more 
distant orbits for captured particles. 

We changed our bulge-disk model by adding a 
smoothed third component - a halo. The structure and the 
shape of dark halos in early-type galaxies appears to be, 
at present, an open question (Sackett et al. 1994, Combes 
& Arnaboldi 1996). For lack of unambiguous knowledge 
about the halo shape, we suppose the halo mass distribu- 
tion to be spherical. As usual, halos are characterized by 
isothermal spheres over some radial interval. For simplic- 
ity, the folowing phenomenological cumulative mass profile 
is used to represent a halo in the present study: 

M h (r) = ^r-Oharctanf^-^ (l-e - ^) , (7) 

where Mho is the total halo mass, r c a cutoff radius, and 
ah a core radius. 

Our choice of halo parameters was determined by 
characteristic values obtained by Reshetnikov & Combes 
(1994a) from the ring kinematics investigation. The core 
radius ah was taken as 9 kpc. The value of the cutoff 
radius r cu t is somewhat arbitary and was taken to be 
equal to 3ah- We choose this value to give a reasonable 
halo mass outside the radius comparable with r per . We 
have let the dark mass to be equal to twice the luminous 
mass (bulge + disk) inside a radius of 15 kpc. The total 
mass of the bulge and disk have been reduced, that is, 
M b + M d = 0.5 = 5 x 10 10 M o , and we let M b /M d = 0.01, 
ab = 0.5 kpc, &d = 3 kpc. 

The geometry of the encounter was the same as in 
previous numerical experiments. Expecting the capture of 
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gaseous matter on distant orbits and formation of a more 
diluted object, we have increased the total number of par- 
ticles (up to 20 000) as well as their size h up to 375 pc to 
keep the gas treatment as a continuous medium. 

The morphology of the new run is shown in Figure |[ 
By the time of t = 6 after the pericenter passage, the cap- 
tured gaseous matter is azimuthally smeared in a annular 
configuration around the galaxy. This structure gradually 
evolves into a closed ring. We have followed its evolution 
up to t — 20. By this time, the ring was completly closed 
but not quite symmetric. According to Rix & Katz (1991) 
for a ring forming as a result of the gaseous satellite dis- 
ruption, the smearing process takes up to seven orbital 
periods - up to 3-4 Gyrs or two times the final time of our 
run. The total mass gathered in the ring-like structure 
does not significally differ from that obtained in haloless 
runs. The ring has a total mass of about 2.2 x 10 9 M Q . 
To explore the sensitivity of this value to the structural 
properties of the donor galaxy, we have constructed an ad- 
ditional model in which the spiral galaxy has its own halo. 
The parameters of this halo were taken the same as for the 
host galaxy. The mass of the captured gas turned out to 
be slightly smaller - 1.4 x 10 9 M Q . (This value is some- 
what small in comparison with that obtained from obser- 
vations. We omit the discussion of this question, because 
this value is obtained for one set of orbital parameters and 
we did not investigate the role of initial conditions, impact 
parameter, and structure of the donor galaxy. It may be 
impossible to understand the formation of very massive 
rings in the frame of the pure accretion scenario and the 
merging event is needed.) 

The most remarkable feature of this annular structure 
is its size. The ring material lies well outside the luminous 
material of the host galaxy. We estimated the diameter 
of the ring of about 30 kpc. This value is typical of the 
objects of the first group of PRGs (see Table 1). The an- 
nular configuration obtained is rather narrow. Its further 
long-lived evolution will lead not only to azimuthal smear- 
ing of the gas but also to flattening of the radial density 
profile, forming a disk-like structure. Two factors promote 
the formation of rings with extended density distributions 
- the viscosity (see Rix & Katz 1991) and the nonspheric- 
ity of the potential (see Fig. 5 in this paper and Fig.l in 
Katz & Rix 1992). 

4. Discussion and conclusions 

We have analyzed the data on global structural properties 
of all known PRGs, and came to a conclusion that there 
is a correlation between the main characteristics of polar 
rings and host galaxies. There exist two different classes 
of PRGs: bulge-dominated galaxies with internal (relative 
to the optical size of a galaxy) rings and disk-dominated 
systems with extended rings. The sample studied is some- 
what small (8-9 objects), but when we combine our prelim- 



inary inference with the results of numerical experiments, 
a plausible picture of PRG properties emerges. 

We have presented the investigation of the ring- 
forming process by studying the stripping of gas-rich 
galaxy outskirts during the encounter with galaxies of dif- 
ferent structures and the subsequent settling of the cap- 
tured material in the potential of a companion galaxy us- 
ing an SPH method. The simulations have demonstrated 
that the ring sizes depend on galaxy potentials. Under 
the same interaction conditions, rings twisting around 
galaxies with a strong concentration of mass to the cen- 
ter (bulge-dominated systems) are less extended (on the 
average) than that around disk- dominated galaxies. More- 
over, we have found that it is impossible to obtain very 
extended rings (with a diameter of about 30 kpc) taking 
into account the distribution of the luminous mass only. 
To obtain such extended rings, the galaxy mass distribu- 
tion should be less concentrated than observed luminosity 
distribution. Therefore the presence of dark halos is re- 
quired to explain the existence of extended rings. 

One can suppose that the observational dichotomy is 
the result of selection. Indeed, interactions between galax- 
ies accompanied by matter accretion occur with partic- 
ipation of various types of galaxies. But we have more 
chances to create long-lived extended rings when the ac- 
creting galaxy is a gas free, disk-like galaxy with extended 
massive dark halo. So such extended polar rings are ex- 
pected to be observed preferably around specific types of 
galaxies. 

The existence of the observational dichotomy suggests 
that PRGs with bulge-dominated central galaxies may 
possess less pronounced dark haloes in comparison with 
the first group of objects (another possibility - bulge- 
dominated galaxies have more centrally-concentrated ha- 
los - docs not change our results as follow from Sect. 3. 1.3). 
These conclusions do not contradict the present-day obser- 
vational data about dark halos in early-type galaxies (e.g., 
Bcrtin & Stiavelli 1993, de Zeew 1995). Nevertheless, we 
do not exclude other explanations of the revealed obser- 
vational dichotomy (for instance, environmental influence 
etc.). 

As was mentioned in Section 2, there is a close simi- 
larity between extended rings and disks of late-type spiral 
galaxies. Extended rings have the same total luminosi- 
ties, integral colors, ratios of HI mass to the blue luminos- 
ity corrected for internal absorption, gas-to-dust ratios as 
normal spiral galaxies; there is evidence for the gaseous 
medium in polar rings to be near (or above) the gravita- 
tional stability limit, as the gaseous disks in spiral galaxies 
(Reshetnikov et al. 1994). Two polar rings (in UGC 7576 
and UGC 9796) demonstrate Ha-derived global star for- 
mation rates to be typical for normal spiral galaxies; gen- 
eral properties of HII regions in these galaxies are nor- 
mal for late- type spirals (Reshetnikov & Combes 1994a). 
HII regions in the ring of prototype polar-ring galaxy 
NGC 2685 demonstrate nearly solar oxygen abundances 
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Fig. 5. Ring formation history for a galaxy with a massive halo. Dimensionless time (counted out from the moment of the 
minimum approach of galaxies) is given in the upper right-hand corner of all frames. The center of each frame coincides with the 
center of the ring-host galaxy. All frames show the orbital plane projection of the ring. The edge length of the frame corresponds 
to 80 kpc. 



(Eskridge & Pogge 1994). Extended rings show radial color 
gradients with bluer colors at larger radii (Reshetnikov et 
al. 1994, Arnaboldi et al. 1995). The HI/H 2 mass ratio in 
the rings of NGC 660 and NGC 4650A is usual for spiral 
galaxies (Combes et al. 1992, Watson et al. 1994). There is 
evidence for the spiral structure in the rings of UGC 7576 
(Reshetnikov & Combes 1994b) and NGC 4650A (Arn- 
aboldi et al. 1996). There is an indication of the pres- 
ence of a large-scale magnetic field in the ring of NGC 660 
(Reshetnikov & Yakovleva 1991). Moreover, according to 
Combes & Arnaboldi (1996) and Arnaboldi et al. (1996), 
dark matter in PRGs could co-exist with the HI com- 
ponent, leaving its mark on the kinematics of the polar 
ring itself similar to the situation with a late-type galaxy 
gaseous disk. One can note finally that PRGs lie close 
to the Tully-Fisher relation for spiral galaxies (Knapp et 
al. 1985). Therefore, PRGs with extended rings could be 
undistinguishable from normal spiral galaxies at less ad- 
vantageous orientation (with a more face-on ring) - see the 
discussion of this question in PRC. Probably, they will 
look like early-type spiral galaxies (maybe barred) with 
extended low surface brightness disks. 

The striking similarity between extended rings and 
spiral disks suggests that PRGs are giant low surface- 
brightness spiral galaxies with decoupled bulges. Although 
this cannot be completely ruled out in some special cases, 
we suppose that the ring capture due to external accre- 
tion is a more frequent process (see Section 1 about ob- 
servations of forming rings). Therefore, in contrast to the 
usual opinion that galactic interactions lead to fast mor- 
phological changes towards SO/E, the observational data 
on PRGs and our modelling simulations provide the ex- 
ample of a rare opposite shift - from an early-type gas-free 
galaxy through the capture of a gaseous disk to a spiral 
galaxy. One can suppose that PRG formation due to ex- 
ternal accretion can imitate in some aspects the formation 
of disk galaxies in the hierarchically clustering model of 
the Universe (e.g., Steinmetz 1996). 

In conclusion, let us reply shortly to the main questions 
of the article (see Section 1): 

the accretion of matter from a gas-rich companion is a 
quite effective mechanism for polar-ring formation; 

the spatial size of the forming ring is determined (un- 
der the same interaction conditions) by the mass distribu- 
tion of the galaxy; 



the process of ring formation takes approximately (7 — 
9) x 10 8 years for haloless systems and reaches up to a 
few Gyrs for host-galaxies possessing massive halos; 

the total mass captured into the ring during an en- 
counter with a gas-rich spiral galaxy exceeds 10 9 M . 
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